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Abstract 

We calculate in (two-loop) chiral perturbation theory the local NN-potentials gener- 
ated by the three-pion exchange diagrams with one insertion from the second order chiral 
effective pion-nucleon Lagrangian proportional to the low-energy constants c\ 2,3,4- The 
resulting isoscalar central potential vanishes identically. In most cases these 37r-exchange 
potentials are larger than the ones generated by the diagrams involving only leading or- 
der vertices due to the large values of 03,4 (which mainly represent virtual A-excitation) . 
A similar feature has been observed for the chiral 2-7r-exchange. We also give suitable 
(double-integral) representations for the spin-spin and tensor potentials generated by the 
leading-order diagrams proportional to g\ involving four nucleon propagators. In these 
cases the Cutkosky rule cannot be used to calculate the spectral-functions in the infinite 
nucleon mass limit since the corresponding mass-spectra start with a non-vanishing value 
at the 37r-threshold. Altogether, one finds that chiral 37r-exchange leads to small correc- 
tions in the region r > 1.4 fm where Itt- and chiral 27r-exchange alone provide a very good 
strong NN-force as shown in a recent analysis of the low-energy pp-scattering data-base. 
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Over the last years effective field theory methods have been successfully applied to the two- 
nucleon system at low and intermediate energies Jl], |2|, ^ |], jl| ]7j . The idea of constructing 
the NN-potential from effective field theory was put forward by Weinberg Jl| and this was 
taken up by van Kolck and collaborators @. Their NN-potential contained one- and two-pion 
exchange graphs based on time-ordered perturbation theory (regularized by a Gaussian cut- 
off) and a host of contact interactions. The method of unitary transformations was used by 
Epelbaoum, Glockle and Meifiner to construct an energy-independent NN-potential from 
the effective chiral pion-nucleon Lagrangian. Based on one- and two-pion exchange and a few 
contact interactions which contribute only to S- and P-waves a good description of the deuteron 
properties as well as the NN phase-shifts and mixing angles below 7] a j = 300 MeV was found 
in that framework ||. 

The issue of the chiral 37r-exchange was first addressed by Pupin and Robilotta || who 
considered one particularly simple diagram with no internal nucleon propagators. In a recent 
work H we have started to calculate the static NN-potentials generated by the (two-loop) 3n- 
exchange diagrams with all possible interaction vertices taken from the leading-order effective 
chiral 7r7V-Lagrangian. In ref.|| we restricted ourselves to the evaluation of four relatively 
simple classes of diagrams with the common prefactor g\/f%. It was stressed in ref. that 
the chiral 37riV./V-contact vertex employed in ref.|| depends on the choice of the interpolating 
pion-field and therefore one has to consider representation invariant classes of diagrams by 
supplementing graphs involving the chiral 47r-vertex. Another obvious consequence of this 



is that there is no unique definition of the so-called "pion-nucleon form factor" which is often 
introduced in phenomenological models. The only meaningful separation of the NN-potential is 
the one into point-like l7r-exchange and 37r-exchange. The four classes of 37r-exchange diagrams 
evaluated in ref. |]] gave rise only to isovector spin-spin and tensor potentials with a strength 
of ±0.1 MeV or less at r = 1.0 fm. Compared to the chiral 27r-exchange NN-potentials || [f) 
these are negligibly small corrections. 

In a subsequent work |L0| the static NN-potentials generated by the 37r-exchange diagrams 
proportional to g\ and g\ (grouped into classes V, VI, VII, VIII and IX) have been worked out. 
The resulting coordinate space potentials have typical strengths of a few MeV at r = 1.0 fm. 
A surprising result of ref.|It| was that the total isoscalar central 37r-exchange NN-potential 



vanishes identically. The static 37r-exchange NN-potentials have been constructed in refs. || pL0| 
from the mass-spectra (i.e. the imaginary parts) of the respective two-loop diagrams in the 
infinite nucleon mass limit M — > oo. This method could however not be successfully applied to 
the spin-spin and tensor potentials generated by the diagrams proportional to g\ (classes VIII 
and IX) involving four nucleon propagators. In these cases the integrand of the 37r-phase space 
integral / § z %<\ duJidijj 2 . . . developed a non-integrable singularity of the form (1 — z 2 ) -3 / 2 on the 
boundary z 2 = 1. By considering the analytically solvable example of the crossed three-boson 
exchange with scalar coupling in ref. |10|| the failure of the spectral-function method could be 



understood from the fact that in the limit M — > oo the mass-spectrum starts with a non- 
vanishing value at the 3-boson threshold \i = 3m. In coordinate space this feature manifests 
itself large-r asymptotics of the form e r r . 

The purpose of this work is to present suitable double-integral representations for these 
remaining 37r-exchange spin-spin and tensor NN-potentials proportional to g\ together with 
numerical results. Furthermore, will work out here the next-to-leading order corrections to the 
chiral 37r-exchange NN-potential. In general, the next-to-leading order contributions are given 
by all diagrams with one insertion from the second order chiral 7riV-Lagrangian C\^. However, 
from an analogous calculation of the chiral 27r-exchange in ref. || one has experienced that the 
dominant contributions come from the second order chiral 7T7riViV-contact vertex proportional 
to the low-energy constants 01,2,3,4, whereas the relativistic 1/M-corrections are much smaller. 
The exception to this is, of course, the spin-orbit NN-potential which is a truly relativistic 
effect vanishing in the infinite nucleon mass limit M — > 00. Therefore we will consider here 
only the two-loop 37r-exchange diagrams with exactly one insertion proportional to 01,2,3,4- The 
corresponding effective chiral 7riV-Lagrangian reads 

C^l = N{2 Cl ml(U + U r ) + C2U0 + c z u v u u + % -c±o ■ (u x u)}n , (1) 

with the axial vector quantity u u = d u £}. The unitary matrix £ 2 = U = l + ir-7r/f w — 
7? 2 /2/ 2 + ... collects the Goldstone pion-fields and f n = 92.4 MeV denotes the weak pion 
decay constant. The low-energy constants Ci.2,3,4 have been determined in a one-loop chiral 
perturbation theory analysis of low-energy 7riV-scattering in the physical region JTTJ] and inside 



the Mandelstam triangle |l2j]. Their central values are c\ = — 0.8GeV 1 , Ci = 3.2 GeV 1 
c 3 = —4.7 GeV -1 and c 4 = 3.4 GeV -1 . We prefer to use the values obtained in ref.|T2" 



since 



at the small sub-threshold pion-energies relevant inside the Mandelstam triangle the chiral 
expansion of the 7riV-scattering amplitude is expected to work best. As outlined in ref.jnj the 
large low-energy constants 02,3,4 receive dominant contributions from virtual A(1232)-isobar 
excitation. In ref. [0] it has been demonstrated that the 27r-exchange NN-potentials generated 
by explicit A(1232)-isobar excitation are very similar to those obtained from the iViW7r-contact 



interaction proportional to 01,2,3,4. In the so-called small scale expansion of ref.|14|| the AN 
mass-splitting A = 293 MeV is treated as a further small expansion parameter and 27r-exchange 
graphs with A-excitation are counted as leading order ones. This is another reason to focus 



on the 37r-exchange potential proportional to 01,2,3,4. The second order chiral 7T7riViV-contact 
vertex following from eq.(l) reads 



y 2 { - 25 ab [2 Cl ml + c 2 q° a q° b + c 3 q a - q b ] + c 4 e abc r c a ■ (q a x q b )} , (2) 

J 7T 

where both pion four-momenta q a and qb are out-going ones. Note that the isoscalar part in 
eq.(2) is spin-independent while the isovector part is spin-dependent. 

Next, let us give some basic definitions in order to fix our notation. In the static limit and 
considering only irreducible diagrams the on-shell NN T-matrix takes the following form 

T NN = V c (q) + V s (q)ai- a 2 + V T {q)& 1 - q a 2 ■ q 

W c (q) + W s (q) ai-&2 + W T (q) a x -qa 2 -q n ■ T 2 , (3) 



+ 



where q = \q\ denotes the momentum transfer between the initial and final-state nucleon. The 
subscripts C, S and T refer to the central, spin-spin and tensor components, each of which 
occurs in an isoscalar (Vc,s,t) and an isovector version (Wc,s,t)- As indicated, the (real) 
NN-amplitudes Vc(q), ■ ■ ■ , Wt(q) depend only on the momentum transfer q in the static limit 
M — > 00. We are here interested only in the coordinate-space potentials generated by certain 
diagrams in which three pions are simultaneously exchanged between both nucleons. For this 
purpose it is sufficient to calculate the imaginary parts of the NN-amplitudes, Im Vc.^t^aO 
and Im Wc,s,T{ifJ>), analytically continued to time-like momentum transfer q = i/j, — + with 
/i > 3m w . These imaginary parts are then the mass-spectra entering a representation of the 
local coordinate-space potentials in the form of a continuous superposition of Yukawa-functions, 

Voir) = --i- r dfifie-^hnVciitJ.), (4) 

~ 1 r°° r 1 

V s (r) = — dfifie-"" fM 2 lmV T (ifx)-31iaV s (iiJ,) , (5) 

V t( t ) = ^~2~s / dfifie-^ r (3 + 3fir + fi 2 r 2 )lmV T {ifi). (6) 

The isoscalar central, spin-spin and tensor potentials, denoted here by Vc,s,t{t), are as usual 
those ones which are accompanied by the operators 1, a\-a 2 and 3 fi\-r o 2 -r— a\ -a 2 , respectively. 
For the isovector potentials Wc,s,t{t) a completely analogous representation holds. The imag- 
inary parts entering eqs. (4,5,6) are calculated from the pertinent 2-loop 37r-exchange diagrams 
as integrals of the iV./V — ► 3tt — > NN transition amplitudes over the Lorentz-invariant 37r-phase 
space making use of the Cutkosky cutting rule as explained in refs. |TI| [HJ . The sub- leading 
37r-exchange considered here occurs at fifth order in small momentum (or chiral) expansion. Of 
course at this order there are also many two-loop 27r-exchange diagrams proportional to ci.2,3,4. 
Since the corresponding spectral-functions will (in general) contain one-loop sub-divergences 
(related to the fourth order chiral 7riV-amplitude) an adequate discussion of these 27r-exchange 
potentials requires good knowledge of the strength of the counter-terms in C^ N jrO). With 



these included the corresponding 27r-exchange potential is well-defined, however, its explicit 
computation goes far beyond the scope of the present paper. 

We start with the diagrams of class X shown in Fig.l. The heavy dot symbolizes the second 
order chiral 7r7r./ViV-contact vertex eq.(2) proportional to 01,2,3,4. Diagrams for which the role of 
both nucleons is interchanged are not shown. These lead to the same contribution to the NN- 
potential, i.e. a factor of 2. As stressed in ref.0 diagrams involving the chiral 3tcNN- vertex or 
the chiral 47r-vertex depend on an arbitrary parameter and therefore one should consider the 
full class X as one entity. Obviously, the last two pion-pole diagrams in Fig.l contribute via 
coupling constant renormalization also to the point-like l7r-exchange. This effect is however 
automatically taken care by working with the physical (on-shell) 7riViV-coupling constant g^N- 



Fig.l: 3n-exchange diagrams of classX proportional to g\. Solid and dashed lines represent 
nucleons and pions, respectively. The heavy dot symbolized an insertion from second order 
chiral 7rN -Lagrangian. The combinatoric factor of these diagrams is 1/2. Diagrams for which 
the role of both nucleons is interchanged are not shown. They lead to the same NN-potential. 

Note that we are considering here only the finite-range Yukawa-parts of the l7r-exchange 
and we disregard all zero-range <5 3 (r*)-terms. Transformed into momentum space the latter 
become polynomials in q 2 with possible contributions from higher- derivative operators. From 
an inspection of the spin- and isospin factors occurring in the diagrams of class X one finds 
immediately that only non-vanishing isovector spin-spin and tensor NN-amplitudes Ws,t will 
be obtained. We find the following imaginary parts from class X (dropping from now on the 
argument i/i), 
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where we have introduced the abbreviations Ci 3 (w) = Ac^ml + c 3 (w 2 — 2ml) anc l ^( w ) = 
w 4 + fi 4 + m\ — 2w 2 fi 2 — 2w 2 ml — 2fi 2 ml . The variable w denotes the invariant mass of a pion- 
pair and its kinematically allowed range is 2m 7T < w < fi — m n . The c?w-integrals in eqs.(7,8) 
could of course be solved easily in terms of square-root and logarithmic functions. However, 
we want to avoid the resulting rather lengthy expressions. 
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Fig. 2: 3n-exchange diagrams of class XI proportional to g\. These give a vanishing contribution 
to the NN-potential. 

Next, we consider the diagrams of class XI shown in Fig. 2. Each of them gives a vanishing 
contribution to the NN T-matrix for the following reason. Since the leading order (Weinberg- 
Tomozawa) 7r7riViV-vertex is of isovector nature (i.e. proportional to e abc r c ) a non-zero isospin- 
factor is obtained only from the spin-dependent isovector term in eq.(2) proportional to c 4 . 



When combined with the spin-independent Weinberg- Tomozawa vertex on the other side of the 
bubble-type sub-diagram one finds immediately that the pertinent 1-loop integral is zero in the 
heavy baryon formalism (basically because the spin- vector a has no time-component). 





* ■ 



Fig. 3: "in-exchange diagrams of class XII proportional to g\. For further notation, see Fig. 1. 

Next, we consider the diagrams of class XII shown in Fig.3. The isoscalar contribution 
comes exclusively from the C4-term in eq.(2). Altogether one obtains the following imaginary 
parts of the isoscalar and isovector spin-spin and tensor NN-amplitudes from class XII, 
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Fig. 4: 37r -exchange diagrams of class XIII proportional to g A . For further notation, see Fig. 1. 

Next, we turn to the diagrams of class XIII shown in Fig.4. The technique to separate off the 
iterative part from the first two reducible diagrams has been explained in ref.0]. The iterative 
(or reducible) part is defined here (entirely) within perturbation theory as that part which 
carries in the numerator the large scale enhancement factor M, the nucleon mass. Such an 



iterative part does therefore not obey the naive chiral power counting rules (see also section 4.3 
in ref.|J on the so-called iterated l7r-exchange). A detailed discussion of possible ambiguities 
showing up in (non-perturbative) iterations to infinite orders (e.g. via a Schrodinger equation) 
can be found in ref . ]T7| . The isovectorial spin-dependent contact vertex proportional to c 4 (and 
in fact only this one) produces now also a central NN-amplitude. Interestingly, its isoscalar 
and isovector components come with a fixed ratio. The corresponding imaginary parts read 
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Furthermore, the isoscalar and isovector spin-spin and tensor NN-amplitudes generated by the 
diagrams of class XIII have the following imaginary parts, 
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Fig. 5: 3tt -exchange diagrams of class XIV proportional to g\. For further notation, see Fig. 1. 



Finally, we consider the (irreducible) diagrams of class XIV shown in Fig. 5. In this case 
the isovectorial C4-term in eq.(2) does not make a contribution to the isovector NN-amplitudes 



Wc,s,t- One obtains from class XIV an isoscalar central NN-amplitude which is however exactly 
canceled by the one coming from class XIII, 



lmVi XIV) = -lmV^ XIII) . (18) 

Therefore, one has the interesting result that there is no isoscalar central NN-potential from 
chiral 37r-exchange even at next-to- leading order (considering the dominant ci,2.3,4-terms). For 
the imaginary parts of the spin-spin and tensor NN-amplitudes generated by the diagrams of 
class XIV one obtains the following expressions, 

Im V^ XIV) = (4/ ^ 2/i 3 J^dw^Jw 2 - Ami - 5/i 6 - 13/i 4 m 2 + 33/i 2 mt - 15m° 

-w\23fi 2 + 37m 2 n ) + w 2 (45m 4 + 8/i 2 m 2 + 19/i 4 ) + 2m 2 (/i 2 - m 2 ) 3 ^ 2 } (19) 
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This completes the presentation of analytical results. In Table 1, we present numerical results for 
the coordinate-space NN-potentials generated by the 37r-exchange graphs of classes X, XII, XIII 
and XIV for inter-nucleon distances 0.8 fm< r < 1.4 fm. We use the parameters f n = 92.4 MeV, 
m n = 138 MeV (average pion mass) and the central value of the nucleon axial-vector coupling- 

H|. For the second order low-energy constants 01,2,3,4 we use the 
3.2 GeV" 1 , c 3 = -4.7 GeV" 1 , c 4 = 3.4 GeV" 1 mentioned above. 



= 1.267 ±0.004 
-0.8 GeV" 1 . 



constant qa 

values Ci = — U.Siorev C2 
One observes from Table 1 that the next-to-leading order 37r-exchange NN-potentials are often 
larger than the leading order ones (comparing e.g. values at r = 1.0 fm). The reason for this 
are the large values of the low-energy constants C3 and C4 representing the strong A-excitation 
effects. A similar feature has also been observed for the chiral 27r-exchange NN-potential in 
ref.||. In that case the effects occured however more selectively in the various channels. For 
example, there was no 27r-exchange isoscalar central potential at leading order and at next-to- 
leading a large isoscalar central potential with a strength of about —300 MeV at r = 1.0 fm was 
generated mainly by the 7T7riViV-contact vertex proportional to C3. Note that the low-energy 
constant C3 is related to the nucleon axial polarizability || The repulsive isovector central 

1.0 fm is also fairly 



potential W^ IH \r) from class XIII with a strength of almost 50 MeV at r 



large. The largest potential is however the attractive isovector tensor potential Wj, 



(XIII), 



At 



distances r < 1.1 fin it actually overwhelms the repulsive isovector tensor potential of the 
(point-like) l7r-exchange. (Since the eigen-values of f\ ■ T2 are —3 and +1 the l7r-exchange 
tensor potential can also act attractively, e.g. in the 3 Si-state relevant for deuteron binding.) 
It is well-known that phenomenology requires a reduction of the l7r-exchange NN-tensor force 
at intermediate distances. The effect obtained here is presumably too strong and in addition 
the isovector tensor potential gets reduced by the chiral 27r-exchange (see ref. [[J). Because 
of the finite size (quark-substructure) of the nucleon one can trust the present calculation 
based on point-like chiral 7riV-interactions at most for distances r > ro — 1 fm. Clearly, as 
soon as the nucleons start to overlap substantially the picture of an undisturbed multi-pion 
exchange between nucleons breaks down. The limiting distance scale r — 4/m p ~ lfm (with 
m p = 770 MeV the p-meson mass) should be considered as a conservative estimate based on 



the argument that the nucleon axial radius 
from chiral pion-loops |I9 



\r6/m p receives no substantial contributions 



Only by confronting chiral multi-pion exchange potentials with 
NN-scattering data one is able to determine this limiting distance scale r more precisely. 
By analyzing the elastic proton-proton scattering data-base below 350 MeV laboratory ki- 



netic energy it has been shown recently in ref . |20| that I71"- plus chiral 27r-exchange gives a very 
good strong NN-force at least as far inwards as r = 1.4 fm (i.e. the length scale corresponding to 
the pion Compton wave-length l/m n ). For such distances (r > 1.4 fm) the chiral 37r-exchange 
NN-potential calculated here and previously in refs.[f|, |l0|) turns out to be a negligibly small 
correction. By inspection of Table 1 one observes that 37r-exchange potentials typically decrease 
by about a factor of 10 each 0.3 fm step outwards. 



Fig. 6: ^-exchange diagrams of class VIII proportional to g\. 

Finally, we like to present results for the leading order spin-spin and tensor NN-potentials 
generated by the 37r-exchange diagrams proportional to g\ (classes VIII and IX shown in 
Figs. 6 and 7) which could not be evaluated with the help of the spectral-function method 
in ref. [10]. The spin- and isospin factors of the two diagrams of class VIII (shown in Fig. 6) 



are such that the irreducible part of the first graph and the second graph give rise only to a 
non-vanishing isoscalar central and isovector spin-spin and tensor potential. For the latter the 
spectral-function method lead to non-integrable boundary singularities in the 37r-phase space 
integral. Therefore an alternative calculational technique is required. 

The coordinate-space potential is generally given as the 3-dimensional Fourier transform 
j d 3 qexp(i q-r) ... of the 2-loop momentum-space amplitude constructed from the propagators 
and vertices of a diagram. This requires in total 11 integrations, / d z qdH\dH2- Usually, one 
performs first the integrations over the loop-energies / i and l 02 using residue calculus and one 
obtains this way the sum of energy denominators of all time-orderings of the diagram. In the 
present case this would lead to a complicated rational function of the three (on-shell) pion- 
energies 1^1,2,3 which could not be further integrated analytically. The strategy is to perform 
a Wick-rotation of both loop-energies, loj — > ihj, j G {1,2}, (of course one has to proof 
that this is indeed allowed) and to evaluate the remaining three 3-dimensional integrations 
analytically. The latter is possible because the original integrand factorizes in the variables l\, 



l 2 and h = q — h + h (see also the appendix in ref.fjnj). The result of this involves among 
other factors (originating from the vertices of the diagram) a product of three Yukawa-functions 



exp(— rymj + l 2 j)/4irr. This way one derives the following double- integral representations for 
the isovector spin-spin and tensor potentials from class VIII, 
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W? IH \r) = rf J 9 ** A - e- 3x (l + x)(16 + 26x + 18x 2 + 6x 3 + x 4 ) 
o(lo7r/^) d r' I 

J-00 7T 2 ClC 2 9(^(2 1 V 

+x(10SQ + 611) + a; 2 (4Sn + 3Q 2 ) + 3a; 3 Qn + a; 4 n 2 )] | , (24) 

where we have used the dimensionless variable x = m n r together with various abbreviations: 

a=\A + Ci 2 , P = ^+Ch 7 = Vi + ^-Cs) 2 , (25) 
£ = a + /3 + 7, g = a ^ + ( 5 7 + 7a; n = a/?7 . (26) 

The asymptotic behavior e ~ 3m ^ r r - 2 Q f Wg^ (r) in eqs. (23,24) indicates that the corresponding 
mass-spectra start with a non- vanishing value at the 37r-threshold /x = 3m n . For the numerical 
evaluation of the principal- value integrals in eqs. (23,24) it is advantageous to employ the identity 
dC/(C)C -1 = Io°° d C [/(C) - /(-OR -1 , where the right hand side is singularity-free. As 
a check we have derived an analogous double-integral representation for the isoscalar central 
potential V^ HI \r) and we have found perfect agreement with the numerical results of ref.fllO 
based on the simpler spectral-function method. 

Lastly, we come to the diagrams of class IX shown in Fig. 7. For the isoscalar spin-spin and 
tensor NN-amplitudes the relevant sum of energy denominators is (uf + w|)/(2u; 4 w|u; 2 ). Since 
its denominator factorizes in the three variables l\, l 2 and I3 = q — l\ + l 2 one can perform the 
corresponding 9 integrations analytically and one gets the following closed form expressions, 

V { s IX) (r) = T^f^I — ( 2 + 2x + **) (2 - 2x - 4a; 2 - x 3 ) , (27) 

V T X \r) = tt/tW — (2 + 2x + x 2 )(2 + x-x 2 - x 3 ) . (28) 
(167T/ 2 ) 3 r 7 

We have corrected here in eq.(28) a small misprint (factor 2) which has occured in eq.(24) of 
ref . [pij . The same features concerning factorization hold for the isovector central NN-amplitude 
and one reproduces the result w£ X \r) = (^/2048vr 3 /> 7 ) e" 3a; (l + x)(4 + 5x + 3x 2 ) of ref. 0. 
For the isovector spin-spin and tensor potentials generated by class IX one repeats (in a slightly 
modified form) the technique used previously for class VIII and one ends up with the following 
representations, 

ux),, _ '/a U„ , , „„ j> [+<°<KsK2 



X 



e - * E (l2ar 2 + l2x~ l Y> + 2£ 2 + 8Q + 12(3 2 + x((2S + 12/?)Q - 6n) 



+x"(2Sn + 2 {a 2 + 3a 7 + 7 2 )) + Ax\Q - a 7 )n + x A Ii 2 



(29) 



W^ X) (r) 



0* 



e -3^( 2 + 2x + x 2 )(2 + x - x 2 - x 3 ) 
e -* s (2E 2 - 4Q + x(2Sg - 6n) 



3(167T/2)3 r 7 

too 71X1(2 d(i 
+x 2 {Q 2 + 2«V) + x 3 {Q + 2a 7 )n + x 4 n 2 



Again, an analogous (double-integral) representation for the isoscalar central potential V c 



(30) 



(ix), 
c l 



turned out to be in perfect numerical agreement with the result of ref . |Hj based on the spectral- 
function method. 



Fig. 7: 3ir-exchange diagrams of class IX proportional to g\. 

Numerical values for the spin-spin and tensor NN-potentials generated by the diagrams of 
classes VIII and IX are presented in Table 1. With a typical strength of a few MeV (or less) 
at r = 1.0 fm they are comparable to the other leading order chiral 37r-exchange potentials 
calculated in ref. . 

In summary, we have completed in this work the calculation of the leading order chiral 
37r-exchange NN-potentials. Furthermore, we investigated the dominant next-to-leading order 
contributions. The (analytical) results presented here and in refs. || [TIJ are parameterfree and 
in a form such that they could be easily implemented in a future empirical analysis of low- 
energy elastic NN-scattering. This way one would be able to determine the maximal range of 
validity r > r of chiral multi-pion exchange in the NN-interaction. 
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